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SUMMARY 



This report deBcribes the "behavior of the heat meter 
when It is used in thermal circuit a which are under the in- 
fluence of transient potentials. The transient heat trans- 
fer systems considered are as follows: 

I. A heat meter suddenly placed upon a hot surface 

II. A heat meter mounted upon the Interior surface of 

the cabin wall of an uninsulated airplane, which 
is climbing through air of decreasing temperature 

III. A heat meter mounted upon the interior surface of 

the composite cabin wall of an Insulated air- 
plane, which is climbing through air of decreas- 
ing temperature 

Analytical and graphical solutions for cases I and II 
are presented. The analytical solutions were derived for 
idealised systems composed of lumped resistances and capac- 
itances: the graphical solutions were effected by the Schmidt 
method (appendix A). In general, good agreement was obtained 
between the analytical and graphical solutions. 

The solution of the differential equations which accu- 
rately describes a composite thermal system of distributed 
resistances and capacitances 1b difficult bo obtain; thuB, 
for case III only the graphical solution is presented. 



UACA ARE Ho. 4107 



2 



This report Is intended to show how analytical and graph- 
ical solutions for the transient state can "be used to inter- 
pret the results obtained when heat meter readings are taken 
during the climb -of -an airplane." The solutions presented may 
not bo directly applicable to a particular problem; however, 
the techniques used should serve as a guide for analysing 
similar problems. 



IUTHODUCTIOff 



In previous reports (references 1 and 2) the heat meter 
has been desoribed and steady state added oircuit .resistance 
and temperature corrections are presented therein. This re- 
port contains an analysis of tho heat meter when it is used 
in transient hoat transfer systems; tho effect of added cir- 
cuit resistance and heat meter temperature corrections will 
not bo considered for the transient state. In some thermal 
systems the temperature potentials are functions of time. 
Knowledge regarding the manner in which the heat meter re- 
sponds to transients is therefore necessary. Analytical and 
graphical solutions for the transient state can be used to 
intorprot the results obtained when hoat moter readings are 
taken during the climb of an airplane. Tho solution of the 
differential equation whioh accurately describes a composite 
thermal system of distributed resistances and capacitances* 
is difficult to obtain. An idealised thermal system consist- 
ing of lumped resistances and capacitances* readily yioldB to 
description through a Boluble differential equation. The 
distributed resistances should be small compared to the ther- 
mal resistance of the remainder of the circuit in order to 
obtain good approximations with thiB method. In the case of 
a heat meter mounted upon a wall, the method consists of 
treating the wall as a lumped capacitor looated at the center 
of the wall Instead of a distributed capacitor. The thermal 
resistance of the wall 1b considered to be divided equally on 
both sides of this central capacitor. Similarly the capac- 
itance of the heat meter is lumped as a single capacitor at 
the center of the meter, and the thermal resistance of the 
heat meter is split equally to both sides of this capacitor. 
The Schmidt graphical method of solving a transient system 
should be used when the Bystem contains a relatively large 
amount of distributed resistance and capacitance or when, the 
temperature potentials or some of the circuit resistances 



•Distributed and lumped resistances and capaoltanoeB are 
discussed in appendix C. 
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vary in somo irregular manner with time; aa under these con- 
ditions the differential equation describing the system "be- 
comes unwieldy. . _ 

This program of research in the Spectr o-Badiometrio 
Laboratory of the Department of Mechanical Engineering of 
the University of California vqb conducted under the sponsor- 
ship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 

The authors wiBh to express their appreciation to Messrs, 
V. P. Berggren and X. H. Morrln for their assistance in pre- 
paring and checking this' report. 

SYMBOLS 

A area through which the heat is "being transferred, . ft a 
a thermal diffuBivlty, 

hr 

C 1 1/2 thermal, capacitance of heat meter, SiS. 

C m thermal capacitance of heat meter, 
C w thermal capacitance of wall, tjy- 
c^ hoat capacity, ^' uU 



p lb °r 

e base of natural logarithms 

f„ unit thermal convective conductance, 3 x 

C hr ft 3 °] 

3tu 

k thermal conductivity. 



hr ft 3 (°P/ft) 
L a significant dimension in equation (8), ft 

n 1 ,n 3 ,ii 3 ,ii 1 constants obtained by simplifying equations - ■ 

(see appendix B) 

Btu 

q. heat flow from airplane interior into heat meter, — — 

■ hr 
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q_ heat flow from heat meter into airplane wall, 4~ 

hr 

q v heat- flow from airplane- wall to outside air, Hi - 

hr 

q,, rate of heat transfer through heat meter, Hi 

hr 

q„ rate of heat transfer through heat meter at steady state, 

Btu. 
hr 

£ 1 1/2 heat meter resistance plus hoat meter-air lnterfaoo, 

or hr 
Btu 

Op 

R 3 air gap resistance plus 1/3 heat meter resistance, 

i TP V t» 

R 4 1/3 heat meter resistance, ■ ■ 



Btu 



R B l/3 heat meter resistance plus heat meter-air interface 

resistance, — — 
Btu 

B g 1/3 heat meter resistance plus contact resistance plus 

1/2 wall resistance, °?_kr 
' Btu 

H 7 1/3 wall resistance plus outside air-wall interface 

resistance, ■ 

Btu 

T absolute temperature of air, °B 

t lf t a the temperatures of the hot and cold Junctions of 
the heat meter thermopile, °r 

t m heat meter temperature at midsection, °T 

t w surface temperature, °T 

Ujn true airplane speed, ft/sec 

a defined "by equation T a = T c - aQ, °T/hr 

V weight density, lb/ft 3 
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A6 tine increment "between steps in graphical method, hr 

Ax distance increment , ft _ 

6. time, hr 



T a amhient air temperature, °TF 



T a air temperature outside plane, °TF 
T e air temperature inside cahin (constant), °7 

Btu 

f r equivalent unit thermal conductance for radiation,- ■ f g 6 j 
f unit thermal conductance due to radiation- and comrec- 

IT?*** o . 

air gap resistance plus 1/2 heat meter- ra-si«tance> - — 

THERMAL STBTHMS AKD SOLUTIONS 
I. A Heat- Hater Suddenly Placed upon a Sot Surface 

A heat meter initially at air t empera-ture-, T lt is 
suddenly placed against a hot surface. The t ime— temperature 
history of the heat meter under these eondtt l.ons . is . desired . 
Tho thermal system consists of a hot surface, an air gap 
(contact rosistanco), tho heat '-meter , and. a unit thermal- con- 
ductance due to convection and radiation (heat meter-air 
interface resistance). Experimental data wero- taken for^the 
above conditions. Initially the. t-emperaturea of the air, heat 
meter, and air-gap vera at 78"° 3", while the .t-empera-tur e of tho 
hot surface waB at 122.8° !. A thermocouple Located. under ,tha» : 
first lamlnatrion of the heat met-er~~was used to obtain a time— 
temperatura history of the heat ..meter (ehjown-in ft&* l). Tha* . 
temperature of the hot surf aca..incr-ea-sed' -slightly a"ft-e-r~~the 
.heat meter 'waa- placed -upon the surface.. At steady -stata the 
unit - thermal conductance -(f c - -+-* f r ) was -experiment ally- date r- 

mined to he 2 Btu/hr ft a °T. The properties, af .JJie- alement a 
in the thermal- system are as indicated .in., table I. Tha prod-?. . 
uct *Ycp for air i s ■ small.. and., can be 'eansi-dared -to be .sero 

in the air gap for -thia system. 
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A . Analytical determination of the time-temperature 
history of a heat meter suddenly -placed on a hot surface . 



Air 

Sap 



S 
Wall 

A 



Air 
T l 



Haat ueter 



The System 



VVjvV 



tm 



i 



V s 

Beat motor 



Cm 
Hi 
H 3 

tm 

*w 

6 

T i 



thermal capacitance of heat meter, - 



Btu 



1/2 heat meter resistance plus surface resistance, 
air gap resistance plus 1/2 heat meter resistance. 



°]P hr 
Btu 

°J hr 
Btu 



average heat meter temperature, 



hot surface temperature, 
time, hr 

ancient air temperature. 



Equating the heat flow from the vail to the meter to the 
rate of heat storage in the meter -plus the heat flow from the 
meter surface to the surroundings (postulating unidirectional 
flow; i.e., heat flow from edges = 0) 



*w ~ *m 



dt 



'm 



d6 



m *m - Ti 



(1). 



or 
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m dQ V E 1 E a / " E a R 1 



(3) 



or 



dt m 



(3) 



where 



v, Kl E a 51 e c n ^e 3 + eJ. 



The solution of the differential equation (3) is (see refer- 
ence 9): 



.. - s ♦ (i . a) 

P ^ P ' 



(4) 



Substituting for P and ft 



<m 



(- + ~) 
H ? 

/ E x + B a \ 

v e x e 3 ; 



T i - 



/Ei + E a \ 

v e x e 3 ;j 



or 



t, 



vhere 



(5) 



is a time constant which is a funotion of 



the heat meter resistance and capacitance and the air gap and 
heat meter-air interface resistance. 

This equation represents the t omperature-tlme history 
for a heat meter suddenly placed on a hot surface of constant 
temperature. The results are tabulated In table II and plot- 
ted together with experimental points In figure 1. The dis- 
crepancy between the experimental and calculated points is 
mainly due to the fact that the hot plate teaperature - .- 
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Increased slightly due to the Increased thermal resistance 
added by the heat meter. 

• B. graphical determination of the t 1 me- temperature 
history of a heat meter suddenly -placed upon a hot surface. 
Figure 2 exhibits the Schmidt solution of the time-temperature 
history of a heat meter suddenly placed upon a hot surface. 
An explanation of the Schmidt technique for solving transient 
heat transfer problems is given in appendix A. She time- 
temperature history at a point inside the heat meter where the 
thermocouple is located is taken from the Schmidt solution, 
figure 2, tabulated in table III, and plotted with the exper- 
imental data in figure 1, The deviation between the experi- 
mental and Schmidt time-temperature histories is seen to be 
negligible. Experimentally it took 2.1 minutes for the tem- 
perature of the heat meter to reach 95 percent of the steady 
Btate heat meter temperature. The time required for the tem- 
perature of the heat meter to roach 95 percent of the steady 
state heat meter temperature as determined by the two methods 
Is tabulated below. 

Experimental result Analytical result Graphical result 
Cnin) (nln) (min) 

2.1 1.4 2.0 



C. Analytical determination of the time-heat flow 
history of a heat meter suddenly placed on a hot surface . 
When the capacitance of the heat rceter Is lumped into two 
capacitors the following equation of heat flow through the 
heat ----- - - - •-- -■ 



meter as a function of time results: 



I 



_Heat_mgte_r 

V\AVV — I — V'WvVA — [ 



I J_ 



1 



At 



The ratio of the heat flow through the heat meter at any 
time 6 to the steady state heat flow (8 - infinity) Is 
given by 
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«= l + 



|_n x - n a (m - n a ) C x B 3 Hj 



e 



+ 



|n x - n a (n a - n x ) C x B 3 E 4 J 




The method need In solving this equation and the defini- 
tion of the terms nj. and n a aro given In appendix B. 

Vhen the values n l§ n 3 , H 3 ( B 4 # B B , O x , and G are 
evaluated for this particular case, the following results: 

-2250 6 

-5.13 e (7) 

Table 17 and figure 3 exhibit the time-hoat flow history 
of the hoat meter as given by equation (7). 

D. Graphical determination of the tlmo-heat flow 
h istory o f a heat meter suddenly nla ced o n a hot surface. - 
Because the Schmidt method yields the temperature throughout 
the heat meter as a function of time, the temperature differ- 
ence across a fixed resistance in the heat meter (as measured 
by the heat meter thermopile) is known as a function of time. 
Tho heat flow-timo history of the heat meter calculated from 
the drop in temperature across the fixed resistance as ob- 
tained by tho Schmidt method is tabulated in tablo V and 
-clotted in figure 3. The Schmidt method reveals that the 
heat flow through the heat meter is within 5 percent of the 
steady state heat transfer rate after 1.4 minutes, while the 
analytical solution yields 1.7 minutes. 



II. Transient Behavior of a Heat Meter Mounted upon the Inside 
of an Uninsulated Cabin Wall of an Airplane In Flight 

An airplane loavea the ground, flying at the rate of 
146 feet per second with respect to the still air, and at an 
angle of 20° with the horizontal. Thus the vertical rate 
of climb is 50 feet per second. The airplane climbs at 
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this rate for 9£ minutes, at which time an elevation of 
37,750 feet has been reached. The temperature of the still 
air on the ground is 60° T, while at an elevation of 27,750 
-feet it- ls--40o j. Some of the atmospheric "data at different 
elevations are tabulated in table VI. (See reference 3.) 
How long must an airplane fly at 27,750 feet, after having 
attained that elevation, "before the rate of heat transfer 
through the heat meter mounted' on the inside of the cabin 
wall is within 5 percent of the steady state rate of heat 
transfer? As the airplane climbs, the air temperature, air 
density, -and the outside unit thermal convectlve conductance, 
vary. The average unit thermal convectlve conductance be- 
tween the air and the outside of the cabin wall may be ex- 
pressed as (of. reference 4): 

The significant dimension chosen in this case is 6 feet. 
(Tub significant dimension L in equation (8) is the length 
of a flat plate over which the average unit thermal convec- 
tive conductance 1b desired. An average unit thermal convec- 
tlve conductance over a 6-foot length was obtained.) By us- 
ing the data tabulated in table VI (reference 3), the unit 
thermal convectlve conductance is calculated as a function of 
elevation and time of flight as phovn in table VII. The 
properties of the elements of the thermal system are as indi- 
cated in table VIII. The unit thermal conductance due to 
convection and radiation on the inside of the cabin is postu- 
lated to be 1 Btu/hr ft 3 °T . 

A. Analytical determination . 



Symbols 

B 3 integration constant 

a constant defined in equation (B?) 

B 4 integration constant 

b constant defined in equation (58) 

c constant defined in equation (BP) 

n 3 ,n 4 constants defined in equation (Ell) 
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first lumped thermal reBi stance, Inside surface resist- 
ance plus 1/2 heat meter resistance, °p/ Btu 

hr 

S fl second lumped thermal resistance, 1/2 heat meter resist- 
ance plus contact resistance plus 1/2 wall resistance, 

*' hr 

t w average wall temperature, °r 

a coefficient of temperature drop with time, — = constant 

hr 

The System 
fWVi 



Cabin air 

T„ 



»1 



Air 
gap 



*2 



1 1 

Heat 
tieter 



J 3 



Outside air 



T 2 = T c " ae 



t w and tjn are 
located at the center 
of the rail and heat 
meter, respectively 



Cabin 
wall 



Thermal circuit 




Heat muter Cabin wall 
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The outside air temperature, T a , varies linearly with time 
and may be expressed as T a a T c . 08. (a Is a constant.) 

-The' complete analytical solution 1b contained in appen- 
dix B. The differential equation obtained Is: 

d6 L >. R 6 e 7 / \ B.i Be 'J d6 

+ *■ * t m + Ba + M t o . qe ( 9 ) 

^ Bx E B B7 ' \ Ex E s B 7 / Be B? 

(ThlB equation and equation (B6) in appendix B are identical.) 
The rate of heat transfer through the heat meter as a func- 
tion cf time is 

^ a ae b a B 2 B 6 B 7 + 0 m a B x 



Bj + B 6 + E 7 (Bi + B e + E 7 )° 2(B X .+ E s + B 7 ) 
n 3 e 

a e T + b n 4 B t B e B 7 _ C m R 1 

KEj + E 6 + B 7 ) L Hi + B B + S 7 2 J 



(n 3 - n 4 

n 4 e 



a 4 0 

ae f b n 3 E, E 
1 + i 

(n 4 - n 3 )(H 1 + H 6 + R v ) L E 2 + B 6 + 



9 il _ °° *i " til ( 10 ) 
R-7 2 J 



(This equation and equation (B26) are identical.) 

Tho physical meaning of this equation can be explained 
as follows: 

The first term represents the heat flow If a steady state 
is reached at each point. The Becond and third terms together 
represent the lag In heat flow reading duo to the capacities 
in the system, after the effect of the Initial conditions has 
vanished. The last two terms represent the effect of the ini- 
tial conditions. The whole equation represents the heat flow 
as read by the meter as a function of. time- 

The numerical evaluation of this equation Is plotted In 
flguie 4 and the calculated values are given In table IX. 
Three curves are present In figure 4. Curve I shows the heat 
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flow through, the wall without the heat meter If a steady 
stflte la reached at each elevation. Curve II gives the heat 
flow reading which would he Indicated by the heat meter if a 
s't'eTady 'atat'e were" reached at each elevation. Curve III gives 
the heat flow ae read by the meter under the assumed condi- 
tions. In figure 8 the calculated values of heat flow are 
plotted with the values ohtalned by the Schmidt method. The 
greatest difference between the two curves is due to the fact 
that an average thermal resistance "between the outside air 
and the wall was employed for the analytical solution; where- 
as, for the graphical solution, the appropriate variable re- 
sistance was used. 

The equations (9) and (10) apply to the climbing air- 
plane, but the heat flow curve after the airplane levels off 
is also plotted. The differential equations and solution 
for this oase are nearly the same as before, the difference 
being that t 8 i b now constant. The initial conditions are 
chosen at the instant the plane levels off. The solution is: 

nj P n4 G , v 

°*M = + B = 6 + B 4 e (11) 

where B 3 and S.i are determined by the heat flow conditions 
at the instant the plane levels off. When the appropriate 
values are substituted, the resultant equation is: 

n = 87.0 - 4.13 e -435 e + 28.83 e ~ 55 * 5 6 (13) 
h 

The calculated values of heat flow are given in table IX end 
are also plotted in figure 4 and represent the part of curve 
III for values of time over 0.155 hour. It is seen that 
af+;er leveling off, it takes approximately 0.035 hour, or 
3.1 minutes for the heat meter reading to reach 5 percent of 
the steady state value. 

Ej Gr aphi cal determination .- The Schmidt solution of the 
t ine-t oraperatur e history throughout a heat meter which was 
attached to the inside of the wall of an eirplane in flight 
was constructed in the following manner. First, the thermal 
resistances of the circuit elerents of the system were repre- 
sented by a convenient aoale. Next, the initial temperature 
distribution was Indicated. Since the thermal resistance of 
the cabin vail was negligible, it was represented by a line 
on the Schmidt solution. The thermal capacitance of the air 
gap was small and neglected. Thus the temperature gradient 
for a particular t ime incremen t was constant across the air 
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gap. The heat meter was divided into one Blah only as the 
temperature gradient^ was rather small. Since Ax = 0.004 foot 
and' a = 0.0 0364 ft /hr, the time increment was 

= 771— — N - 0.00303 hour or 0.18 minute. The outside 
2a 3(0.00264) 

air temperature and the value of k/f c are plotted against 
time in figure 6. .Thus for any- time (determined hy the number 
of time increments) the directional point R can be detert-. 
mined. Table X and figure 8 exhibit the time-heat flow re- 
sults through the heat meter obtained by the Schmidt method 
shown in figure 7. It would take 1.25 minutes for rate of 
heat transfer through the heat meter to reach 5 percent of 
the steady state value after the airplane leveled off. 

The steady state rate of heat transfer is that rate 
which would occur at each particular elevation if the plane 
were to level off and fly at that elevation for an infinite 
time. Prom figure 8 it can be seen that the transient rate 
of l.eat transfer as measured by the heat meter is greater 
thar, the steady state rate of heat transfer. This observa- 
tion can be explained as follows: 

A greater amount of the heat stored in that section of 
the heat meter nearest the decreasing outside air will flow 
ov.t of ' the h.->at motor than that stored in the section of 
the heat meter which is farthest from the outside air. Thus 
the tpraporature of the first section will be decreasing at a 
creator rate than thit of the s.:cor.d section, and since the 
heat meter' thermopile erasures the temperature difference 
between these two sections within the heat motor, the tran- 
sient rnto of heat transfer as measured by the heat meter 
will be greater than the steady state value. 



III. The Transient Behavior of a Eeat Meter Mounted upon the 

InBide of an Insulated Cabin Vail of an Airplane in Flight 

If a heat meter is attached to an insulated cabin wall 
of an air-plane flying under the previously specified condi- 
tions (illustrative example II), how long must the airplane 
fly at 27,750 feet in order that the heat flow through the 
heat meter ae within 5 percent of the steady state heat flow 
when the heat meter is mounted between the wall and insula- 
tion, and x*hen the heat meter is mounted at the inBulat ion- 
cabin air interface? The thermal circuit consists of an out 
side unit thermal conductance, a cabin wall » a heat meter, a 
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contact resistance, an insulation material, and an inside 
unit thermal conductance as shown in figure 9. The distrib- 
uted resistance and capacitance of the heat meter is negligi- 
ble compared to the large distributed resistance and capaci- 
tance of the insulating material. Thus the heat meter very 
nearly measures a rate of heat transfer that would occur if 
the heat meter were not in the thermal circuit . Since the 
system contains a relatively largo distributed resistance and 
capacitance of the insulating material, the analytioal method 
of lumping the circuit constants becomes cumbersome, and the 
Schmidt method of solution is used. The properties of the 
insulating material are: 

k = 0.025 3tu/hr ft 8 — 

ft 

V - 5 lb/ft 3 
Thickness => 0.5 inch 
0 p = 0.40 Btu/lb °T 

k 0.025 _ „_ a 

a = — — = r- = 0.0125 ft /hr 

Cp'V 0.40X.5 ' 

Thermal Circuit 



Outside air- Cabin air- 




heat meters 
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The Sohmidt solution (fig. 10) was carried out as in the 
previous case. The insulation was divided into eight small 
slabs; Ax = 0.0052 foot and the time Increment was 

A0 = (0-0052) = 0.001085 hour or 0.065 minute. Tables 11 
2(0.0125) 

and 12 and figure 11 exhibit the time-heat flow history of 
the heat meter mounted on the cabin wall as denoted above. 
The tiaes required for the heat flow into the insulating ma- 
terial as measured by the heat meter from the cabin and the 
heat flow out of the cabin wall into the outside air to be 
within 5 percent of the steady Btate heat flow value after 
the plane has leveled off, are tabulated below: 



Heat flow into insulation Heat flow out of wall 
(min) (min) 

3.9 4.2 



00IT0LUSIOHS 



1. The analytical and graphical solutions of the partic- 
ular transient heat transfer problems presented in this re- 
port may be used as a guide in analyzing transient heat trans- 
fer problems in general. 

2. If a heat meter is used in a thermal system where any 
of the quantities such as t 3i t 0 , or f 0 vary in some un- 
known manner with time so that it is impossible to obtain 
analytical or graphical solutions, the heat transfer rates as 
measured by the heat meter cannot be readily Interpreted. 

3. In the oase of a heat meter suddenly applied to a hot 
surface of constant and uniform temperature, a heat meter 
reading cannot be made until steady state obtains (when the 
temperature 1b constant with time). The time required for 
steady state to occur can be predicted analytically or 
graphically. 

4. In the case of a heat meter mounted upon the cabin 
wall of an airplane and the temperature of the outside air 
decreasing linearly with time, analytical and graphical solu- 
tions reveal that large errors in determining rates of heat 
transfer can be made by taking readings during the transient , 
and considering them as steady state. For the specific caee 
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presented in this report, the difference between the Bteady 
state and transient heat transfer is directly proportional 
.to the rate., of ohange (a) of the ambient air temperature and 
1 s given "by: 



5. Tor an airplane cabin with an appreciable amount of 
Insulation the thermal circuit is difficult to solve by ana- 
lytical techniques presented because the system contains rel- 
atively large distributed resistances and capacitances which 
do not lend themselves to the analytical method of lumping 
circuit constants. Tor those cases the Schmidt method can 
be used. A knowledge of the thermal capacitances of the in- 
sulation may be obtained experimentally (in flight) by using 
heat meters in the thermal circuit; one could be placed on 
the side of the insulation contiguous to the cabin air, and 
the other could be placed between the Insulation and the 
cabin wall. The amount of heat stored in the insulation may 
be calculated from the difference in the two heat meter read- 
ings . 




University of California, 

Berkeley, Calif., May 1944. 
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APPEND IX A 

TEE SCHMIDT (GRAPHICAL) METHOD OF SOLUTION OP THE 
CONDUCTION EQUATION POE THANSIENT UNIDIHECT I ONAL HEAT PLOW 



Heat conduction problems that are very difficult to solve 
analytically may "be solved graphically by the Schmidt method. 
The differential equation for unidirectional heat flow in a 
solid is 

„ a aft 
ae a* 3 

where t is the temperature at any point, 9 the time, z 
the distance from one face, and a the invariable thermal 
diffusivity of the slab. The graphical integration of the 
above partial differential equation is accomplished by replac- 
ing the differential quantities by small finite differences. 
(See references 5, 6, and 7.) Let A9 represent a small but 
finite increment of time 6, and let Ax represent a small 
but finite increment of distance x. The distance increments 
will be denoted by the subscripts (n-i), n, (n+i), and so 
forth, and the time increments will be denoted by the sub- 
scripts (m-i), m, (m+i), and so forth. Thus t n>m will re- 
fer to the temperature of the slab at the distance nAx from 
the surface at the end of mAG time units. Thus the differ- 
ential quantities are represented in finite difference form 
as follows: 

At = t n, lm+i ) " t n , m 
36 AG 

The temperature gradient after mA6 time units between the 
n t3a and (n+l) th increment is, 

{ b&\ _ t( n +i) ,m ~ t n|B1 
\Ax/ + Ax 

The temuerature gradient after mA6 time units between the 
(n-i)t 11 " and n*^ increment is, 
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The rate of change of the temperature gradient at the center 
of the n* h layer is, 

A 3 * = *(n+i).m + tfn-O.m ~ 8t n.m 

When the finite difference quantities are substituted into 
the oonduotion equation, the following equation results: 

tn, (m+Q ~ t n,m ^ t Xn+i),m + a " ^n.m'N 

Ae " a V Ai a ' 

or 

*».c-i) - »».. - ^ [''"■' ■ »;''»'■■ -'.,.] 

Ax u 2 J 

If Ax and A6 are so chosen that 



3a A 6 , 

"S 5- ■ 



t 



f(m+l) = li*±llt* + liSzlh 



n 



2 



The above equation states that the temperature of any distance 
increment at any time is equal to the arithmetic mean of the 
two distance increments on either side at the previous A@ . 

The distance Increment is chosen on the "basis of the in- 
itial temperature distribution. Tor sudden temperature 
changes the distance increments should he chosen small enough 
so that the broken lines connecting the points which repre- 
sent the temperatures of the distance increments will not dif- 
fer too greatly from the continuous temperature distribution 
which would result if distance increments become infinitely 
small. The corresponding time increment that satisfies the 

equation, 2a A a 9 =1, is used. 
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If the thermal system contains an air-Blah Interface 
resistance, a heat rate balance across the boundary must be 
made so that the temperature at the surface of the slab can 
be determined. 



Thus the slope of the temperature distribution at the bound- 
ary is equal to the difference in temperature between the 
slab surface and the air temperature divided by the ratio 
k/f . The ratio k/f may vary depending upon the variation 
of f . 

The sudden oooling of. an infinite slab with an air-slab 
Interface resistance will be used to illustrate the Schmidt 
construction. The temperature of the air T a , the thermal 

conductivity of the slab k, the thermal diffusivity of the 
slab a, and the unit thermal conductance due to convection 
and radiation f are constant valuea. The initial tempera- 
ture of the slab Is t x . In figure 13 the slab is divided 

into distance increments Ax. If the directrix E, located 
at ordinate T a and at a distance of k/f from the surface 

of the slab, and the point 0 (initial temperature t^ at 

the aurfaco at e = 0) are oonnected by a line HO, the 
temperature gradient at the surface of the slab is equal to 
the slope of the line HO. The intersection of lino £0 

with a plane at a distance ^ from the surface looates a 

point which may be considered as being on the temperature 
curve. The following step consists of connecting points a 
and 3 from which l 1 is obtained. Next, Rl 1 is drawn, 
and 3' Is located by connecting points 1 and 3. Thus 
the temperature distribution after one A0 , 0 1 , 1', 3', 3 1 , 
4 1 , and so forth, has been obtainod. This procedure 1b 
continued for each new A8 . 



or 




f 
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APPENDIX B 

THE SOLUTION OF THE TRANSIENT BEHAVIOR OP A HEAT METER 
MOUNTED UPON THE CABIN OP AN AIRPLANE IN PLIGHT 

The cabin air is maintained at a constant temperature 
Tit and the heat meter, cabin vail, and outside air are ini- 
tially at this temperature'. It 1b postulated that the temper- 
ature of the outside air drops linearly with time. By treat- 
ing the vail and meter as lumped capacitors and lumping the 
thermal resistance as described on page 11, it is possible to 
calculate the temperatures of the heat meter and wall, as 
functions of time. It is alBO possible to calculate the heat 
flow from the cabin into the meter, from the meter to the wall, 
and from the wall to the outside air. Por this case the heat 
flow indicated by the meter is approximately equal to the 
arithmetical average of the heat flow as calculated into and 

q + q 

out of the motor, or q m » 21 a. 

2 

The solution using lumped parametors 1b as follows; 
A heat balance on the moter yields: 

. - IsJ^s „ o m ils + luSJ* (B1 ) 

1 R t d9 E 6 

While a heat balance on the wall yields: 

t - t dt t - t 
n «= -JS St = 0 „ + ~S 2 (B2) 



8 R a * d9 R 7 



where 



t w - T a 

Adding equations (31) and (BS) , gives: 
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From equation (Bl) is obtained: 

Differentiating (54) with respect to 8, gives: 
Substituting (B4) and (B5) Into (B3) giveB upon rearrange- 



ment : 

C G 
m v d 6 s 



+ (li * B fl + M t « R g + E ? + -la— 

V H, B fl R 7 / * B, a fl E ? B 6 H ? 

m /B a j- B 6 + E 7 \ _ _a_e_ (3e) 
V B i H s R 7 / E 6 E ? 

This Is tho differential equation for the heat meter 
temperature as a function of time. To simplify thlB equation, 
let 

a = C m C w <E7) 

D a °° (B « + V + C « < B 1 » H »>^(B8) 
R 6 E 7 H x E 6 

rB- + B* + E 



( 6 7 ) (B9) 

V B i B « E ? / 



Then (B6) is rewritten as follows: 



3 



d6 d 0 \ E x E a E 7 / 



T C - -=J«- (BIO) 



i *e *? ' B o B 7 
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And let 



-d ,+ i/b a 



n 3 . _____^2L__- 4ao 
3a 

and 



B -T> -Vd* - 4ao (BllO 
3a 



The solution of the reduced equation la 

n a 6 „ n 4 6 



t m = A x e"* ° + e" 4 * (Big) 



where Aj. and B x are integration constants dependent upon 
the initial conditions. 

The particular solution is 

t m „ t 0 + LJLjMjb »« a r) (B18) 

E 1 + H fl + E ? (fi 1 + H 6 + B. 7 ) a 
and the complete solution is 

3 +Bfl*+T a B l g 



t_ = A, e + B e + t - 

mi l c 



B ! + B a + P-i 



+ p Q» ^i(B-i B a B7) 

(b x + a a + a 7 ) a 



(B14) 



The initial conditions for determining A x and B 1 

are that when 6 - 0: t m = T, ; QIe. = 0. From eouation 
(BID 46 

Aj + Bj_ = * ft .'l ( »l'oV (B15) 

(a x + B a + a 7 ) 

And from the derivative of equation (314): 
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n 3 A x + n 4 B 1 s a Bl (B16) 



E x + E 6 + E 7 



Solving equations (B15) and (316) for A x and B x gives: 

5 &i [i + I V a i H a V I 

- n 4 )(E a + B 6 + E 7 ) L (E a + E fl + E ? ) a J 



and 

B 



9L*a fi + j n s R i *fl B T ] 

(n 4 - BjJfHj + E 8 + B 7 ) L (E x + E 6 + E 7 ) 3 J 



And on substitution Into equation (311) the equation for 
the temperature variation with time of the heat meter 1b ob- 
tained : 



*m = T 0 - 



tt B-! 6 + b a B 1 (E 1 B fl _B g _ ) 
H a + Eg + E 7 (Ej + E B + E 7 ) a 



n 3 6 



(n 3 - n* )(E 3 + E 6 + E 7 ) L R x + R 8 + Ey J 



+ * ** *** 9 \l * » ** ** ^ 1 (317) 

(n 4 - n 3 )(E 1 + E 6 + E 7 ) L E x + E e + B 7 ' J 

To obtain the equation for the temperature variation of 
the wall of the airplane a similar method is employed. The 
differential equation is nearly the same bb before and is 

+ b *** + ct w - + . B ° + M t c _ * V»e _ (3ib) 

d6 a d6 V H x Eg E 7 / E x Eg E ? B ? 

where a f b, and c are the same as In equation (B6) . 
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The solution of the equition is then 



t = A, 



W ° 3 0 H X + 3 6 + S 7 



s + s e + a, (a a + s 6 + s 7 ) 



a 



(B3.9) 



The initial conditions are, when B ~ 0; t v = T c ; — * = 0 
bo that from equation (BIG) ^ 



A. ♦ B a . -°-3 « »* »*. . * * »«>** »; *' (330) 

Sj + E.g + S7 (-^1 Eg 



and 



n 3 A g + n 4 B g = (. H l * H A )a (331) 



This gives the equation for t ae: 



(S, + 2 6 )(j. e 0 n a E i H e b + H 6^ H l H 6 B 7> 



t = t - - _ _._±i ; 2_ + 

w 1 c T 



P. 5 + E e + S 7 Ej. + E s + S 7 (Ei + E a + E 7 ) 

_ a (l a 4- H a ) e na 6 r + h n, B x B fl E ? _ 0 m » t a n 
(n 3 - n 4 )(Sj + E s + E 7 ) L Ei + E a + E 7 S x + 2 e J 



n 8 

cc. (E + E ) e * r b n B S E n , 0 E E 

+ 1 2_ h + 3167 - 3 m 1 8 ] (B33) 

(n„ - n 1 )(S 1 + S 0 + E ? ) L E 1 + E 6 + E ? E x + E e J 
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Thus far the equations for t m , the heat meter tcr.per- 

Jiture, and t v the cabin wall temperature hare "been obt?i:"*i 

To calculate the heat flov as read "by the heat meter it i.a 
noceesar7 to obtain q x and q a , or the Jjeat flow into and 

out of the meter. From equations (Bl) ani 'BS?): 

q 3 T ° : *» (bi) 



arid 



*m ~ *w 



fl e 

Trom equations (Bl) and (317) 

a 6 15 a R ! H 0 E„ as 



n 3 e 



1 H x + S s + E 7 (R x + H fl + H-) 3 (n 3 - n 4 )(S 1 + S 8 + E 7 ) 

n 4 e 



h n„ E, R„ E 
+ — 
E 



a + E e + S 7 J (n 4 - n 3 )(E x + E a + E 7 ) L E x + E B + E 7 J 

(B23) 

n"d -'i-ov. eauations (32) and (314) and (319): 
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*. . H a (*x + He + By) E a (H a + B a + B 7 ) Eg (Ej^ + E 6 + ^7 ) 

_ T> o(Bj + Bh )B l ^S, + 0 m g R, + o B, 6°° 6 

a„ (Ei + Bg + B, ) a Bi + Bg + By H 6 (na - lO-fRi + He+B, ) 



xu e 



frii4 R-i B 6 Hy "I _ a (Ex + B 8 ) e 3 0 "T Ri B c B7 "1 

E t + B 6 + ByJ E 6 (n 3 - n4)(Bi + E a + Rt)'-L Bx + R 8 + R7-I 

+ a ° m At 8 B i + * B i 8 

(n 3 - n 4 )(B 1 + B e + By ) B a (n 4 - n 3 )(B x + ^ + B ? ) 

n 4 e 



[*i + T>n, Hi P. w R ? I g (R, + R B ) 

L B x + R B + B 7 J R E (n 4 - n 3 )CBj + 



e 

R 6 + R 7 ) 



[1 + bn » Ri a s R ? _ £a c m H i e r 1 (B34) 
L R a + R Q + R 7 Bj, + B a J 

which reduces to: 

g 6 T> g B a B a B ? C m a B^ 

18 " E l + H 6 + *7 " <*1 + *8 + E 7 )S + H l + H 8 + H 7 

!L5l1! fi + » n « H ; »o a 7 - o B b, nj 

(n 3 - n 4 )(Bj + B 6 + B ? ) L B t + B fl + B/ J 

g e n * ® f* "bn B B B 1 , 
_ 1 + ° n a *l a fl a 7 - C m B t n 3 1 (B25) 

<n 4 - n 3 )(B 1 + B 6 + B ? ) L B t + B 0 + B ? J 



and nlnce 
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4, 



a 6 



"b re Bj B e B 7 



O m a, B a 



E x + E 8 + B 7 (B x + B e + B 7 ) 3(Ba + B e + E 7 ) 



e aa 6 r + *n 4 B T B a B 7 _ O n B, n 4 "| 

B x + B 6 + B 7 ) L Ex + B 6 + B 7 3 J 



(n 3 - n 4 )(B 



(n 4 — HaM*! + H a + B 7 ) 



[ 



x + B l H 6 B 7 

B, + E- + E, 



°m *i n 3 

-3 



] 



(B36) 
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APPEHDIX 0 



- - .A d istributed reel atance and capacitanoe in a thermal 
system can "be idealized so that the distributed reBistanoe 
may he oomhined into one or more resistors, and the distrib- 
uted capacitance may he combined into one or more capacitors. 
This technique of combining the distributed resistance and 
capacitanoe will yield a system containing lumped circuit 
par ame t er s . — — 



wviwv 



T 



(a) 



T 



i 



(b) 



(c) 



Distributed resistance 
and capacitance 



Lumped resistances and 
capacitances 



Tor instance, diagrams b and c are examples of elements which 
hare been obtained by lumping the distributed resistance and 
capacitance of element a. The degree to which diagram a 
would be simplified (b or o) would depend upon the elements 
of the complete system and upon the accuracy desired. 
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FROPERTBB OP BLEMBMTS IB TUB TBBBMAL 8TSTBI FOR CASS I 





Theraal 

Conductirlty, k» 
Btu/hr t&i*r/et) 


HMtt 

Capacity, o_ 
Btu/lb * 


Height 

DenaltyF 

lbi/rt 8 


• Thiokneea.ft 


Air 
0»P 


0.016 




0.000766 


Heat 
Meter 


0.18 


0.56 




0.004 



(f Q + t f ) - 2.0 Btu/hr ft 8 °F at steady tftate 
TABU II 



Beaulta of Analytioal Solution 
(Eq. 6. Fig. 1) 



Tim 
minutes 







Eaat Hater 
Temperature 










0 






78 


0.86 






96.6 


0.60 






106.7 


0.76 






111.6 


1.00 






116.1 


1.50 






118.8 


2.00 






119.4 


OO 






119.9 



TABIX III 



Beeulta of Qraphloal 8otaddt Solution 



TlM 


Beat Meter 

Temperature 
Of 


0 


78 


.04 


80 


.10 


86.6 


.22 


94 


.55 


100.4 


.56 


108.4 


.91 


114.2 


1.57 


117 


1.60 


117.6 


1.94 


118 


2.17 


118.2 


OO 


119.9 



L 
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Analytical Time-Heat Flow History of a Heat Meter 
Placed upon a Hot Surface, Soe Fig. 5 



(Equation 7, Page 11) 



Time 
Min. 








Btu/ir ft 2 


0 


0 


0 


0.01 


1.50 


124 


0.02 


2.50 


207 


0.03 


3.16 


262 


0.05 


3.84 


318 


0.07 


4.08 


338 


0.10 


4.08 


338 


0.15 


3.80 


315 


0.20 


3.48 


288 


0.25 


3.12 


258 


0.30 


2.92 


242 


0.50 


2.15 


178 


0.70 


1.68 


139 


1.00 


1.32 


109 


1.50 


1.09 


90 



TABLE V 



Graphical Results by the Schmidt Method 
of the Heat Flow-Time History of 
a Heat Meter Suddenly Placed upon a Hot Surface 

(see Figs. 2 and 3) 



Time 
Min. 


+> 0 


s*- 

°F 




Btu/Kr ft 2 


0 


78.0 


78.0 


0.0 


0 


.0072 


78.4 


78.0 


0.4 


39.0 


.0144 


79.5 


78.0 


1.5 


146.0 


.0216 


80.8 


78.4 


2.4 


234.0 


.0288 


81.9 


79.0 


2.9 


283.0 


.0360 


82.8 


79.6 


3.2 


312.0 


.0432 


83.7 


80.2 


3.5 


342.0 


.0720 


86.6 


82.9 


3.7 , 


36,1.0 


.1008 


89.0 


85.6 


3. "4 


332.0 


.1296 


91.1 


87.9 


3.2 


312.0 


.1584 


93.2 


90.1 


3.1 


302.0 


.1876 


95.0 


92.1 


2.9 


283.0 


.2164 


96.7 


94.0 


2.7 


264.0 


.2452 


98.4 


95.8 


2.6 


254.0 


.2740 


99.9 


97.4 


2.5 


244.0 


.3028 


101.1 


98.8 


2.3 


224.0 


.3316 


102.5 


100.3 


2.2 


214.0 


.447 


106.6 


104.8 


1.8 


176.0 


.562 


110.0 


108.4 


1.6 


156.0 


.677 


112.3 


111.0 


1.3 


127.0 


.792 


114.0 


112.9 


1.1 


107.0 


.907 


115.2 


114.1 


1.1 


107.0 


1.022 


116.1 


115.1 


1.0 


97.6 


1.366 


117.8 


116.9 


0.9 


87.8 


«<» 


120.75 


119.9 


0.85 


82.8 
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Proporties of the Atmosphere 



EleV. 

ft. 


Time of Flight 
minutes 


Air 
Temp. 
°F 


Air Denslty.T 

lb/ft 3 


Air Press. 
lb/in 2 abs. 


0 


0 


60 


0.0765 


14.7 


3000 


1.00 


40.3 


0.0700 


13.17 


5000 


1.67 


41.2 


0.0659 


12.23 


10,000 


3.33 


23.3 


0.0565 


10.11 


12,000 


4.00 


16.2 


0.0530 


9.35 


15,000 


5.00 


5.61 


0.0481 


8.29 


20,000 


6.66 


-12.3 


0.0408 


6.75 


22,000 


7.33 


-19.5 


0.0281 


6.20 


25,000 


8.33 


-30.2 


0.0343 


5.46 


30.000 


10.00 


-48.0 


0.0286 


4.36 


35,000 


11.67 


-65.8 


0.0237 


3.46 


40,000 


13.32 


-67.0 


0.0234 


2.72 



TABLE VII 



The evaluation of f„ as a funotion of altitude and time 



eleva- 
tion 
feet 


time 
nin. 


Tz 

°F 


T 
°H 


Ha 

ft/seo 


Y<Ur 

lb/ft 3 


0.80 
faao* air} 

\Lfl.2sy 


T 0.3 


Btu^ ^ 


k 

?e 


0 


0 


60 


620 


146 


.0765 


4.82 


6.68 


20.0 


0.00*1 


3000 


1 


46.5 


508.4 


146 


.0702 


4.60 


6.60 


is.es 


0.00443 


6000 


2 


37.6 


497.6 


146 


.0641 


■ 4.19 


6.46 


17.6 


0.00471 


12000 


4 


16 


476 


146 


.0532 


3.60 


6.36 


14.6 


0.00666 


18000 


6 


-5.5 


464.6 


146 


.0436 


3.06 


6.30 


12.26 


0.00673 


24000 


8 


-27.0 


433 


146 


.0360 


2.58 


6.80 


10.2 


0.00816 


50000 


10 


-46.0 


412 


146 


.0264 


2.18 


6.10 


a. 49 


0.00971 



f B ■ 0.64 T°-3 Ah? \ °- 80 , quatl6n (,) 

f r on outside Is neglected compared to f c I ifi.ZA I 

TABLE VIII 
PROPERTIES OP THE SYSTEM FOR CASE II 



element 


k op 
Btu/ly ft 2 f£ 


thickness 
ft 


lb/ft 3 


o p 

Btu/^F lb 


heat meter 


0.0824 


0.004 


(1.43)(62.4) - 89.3 


0.3S 


air gap 


0.016 


0.000 533 




0 


aluminum 
cabin wall 


116 


0.00267 


165 


0.21 
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Analytical Results of tha H.«t Flow-Time 
History of a Heat Meter Nountad upon the Inside of the Cabin Wall of an Airplane in Flight 

(Equation 34, Figure 3) 







© 
hours 


e 

ninutaa 


4M Btu _ 


hour* 


e 

minutes 


Th, Btu 

a, nr l & 


0 


0 


0 




0.155 


9.3 


111.7 


.005 


0.3 


5.56 


■ 


.157 


9.42 


111.1 


.01 


0.6 


13.3 




.160 


9.6 


108. 3_ 


.015 


0.9 


20.6 




.162 


9.72 


106.3 


.02 


1.2 


26.2 




.165 


9.9 


103.5 


.03 


1.8 


36.1 




.170 


10.2 


99.5 


.04 


2.4 


44.0 




.175 


10.5 


96.6 


.05 


3.0 


61.0 




.185 


11.1 


92.4 


.076 


4.6 


66.5 




.195 


11.7 


90.1 


.10 


6.0 


81.0 




.205 


12.3 


88. S 


.126 


7.6 


95.0 




.235 


14.1 


87.45 


.16 


9.0 


109.0 


1 .255 


15.3 


87.09 



tabus y 

The Heat Flow-Time History of a Haat Meter 
Mounted upon the Inside of the Cabin Wall of an Airplane in Plight 
by the Sehmidt Method 



& 


*2 


*1 


*2 " *1 




hra 


°F 


°p 


°r 


Btu/hr ft 2 


0 


60.0 


60.0 


0 


0 
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WACA ARR «C *I07 TABUS XI 

Results by the Schmidt Method of Heat Transfer out of tha Coxpoiit* Cabin «kll 
Into the Outside Air at Measured by tha Heat Katar 
A« a Function of Tina 
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TABLE XII 



Results by tha Schmidt Method of Haat Transfer into tha Composite Cabin Ma. 11 
From the Cabin a* Measured by the Heat Meter 

As a Function of Time 
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Figure 1.- Time-temperature history of heat meter 
suddenly placed on a hot surface. 
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Figure 4.- The analytical solution of the time-heat flow history of a heat meter mounted 
on the inBide of the cabin wall of an airplane in flight. 
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Figure 5— The outside thermal conductance due to forced convection on the cabin * 
wall as a function of time. 1? 
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Figure 7.- The Schmidt solution for a heat meter mounted on the inside of an uninsulated cabin wall of an airplane in flight. 
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Figure 8.- The tiaie-heat flow history of a heat meter mounted on the inside of the cabin 
wall of an airplane in flight. 



NACA ABE No. 41.07 



Pig. 9 



t varies linearly 
with time from 
60 Off to -40 o? 
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Figure 9.- The thermal circuit of the heat meter placed upon an 
insulated (composite) cahin wall, 
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Figure 11.- Schmidt solution of the time-heat flow history of a heat meter mounted on the com- 
posite cabin wall of an airplane in flight. 
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Figure 12.- The sudden cooling of an infinite slat with an air-slab inter- 
face resistance. 



